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Abstract 



We analyze tests of electron flavor conservation that can be performed at 
the Sudbury Neutrino Observatory (SNO). These tests, which utilize 8 B solar 
neutrinos interacting with deuterium, measure: 1) the shape of the recoil elec- 
tron spectrum in charged-current (CC) interactions (the CC spectrum shape); 
and 2) the ratio of the number of charged current to neutral current (NC) 
events (the CC/NC ratio). We determine standard model predictions for the 
CC spectral shape and for the CC/NC ratio, together with realistic estimates 
of their errors and the correlations between errors. We consider systematic 
uncertainties in the standard neutrino spectrum and in the charged-current 
and neutral current cross-sections, the SNO energy resolution and absolute 
energy scale, and the SNO detection efficiencies. Assuming that either matter- 
enhanced or vacuum neutrino oscillations solve the solar neutrino problems, 
we calculate the confidence levels with which electron flavor non-conservation 
can be detected using either the CC spectrum shape or the CC/NC ratio, or 
both. If the SNO detector works as expected, the neutrino oscillation solu- 
tions that best-fit the results of the four operating solar neutrino experiments 
can be distinguished unambiguously from the standard predictions of electron 
flavor conservation. 

PACS number(s): 11.30.Fs, 26.65.+t, 13.15. +g, 12.60 -i 
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I. INTRODUCTION 



We assess quantitatively the possibility of detecting electron flavor non-conservation 
using 8 B solar neutrino interactions in deuterium at the Sudbury Neutrino Observatory 
(SNO) JT|. The separate conservation of the lepton (electron, muon, tau) flavors is a well 
known ingredient of the standard electroweak model and of some of its extensions. 

Figures ff] and || (and Table [j] and Table |T|) summarize the power of the SNO experiment 
to find new physics. We urge the reader to look at these two summary figures (and tables) 
before descending into the necessary details, which are analyzed in this paper. 

Solar neutrinos offer a unique possibility to detect electron-flavor non-conserving pro- 
cesses. In solar neutrino experiments, a pure beam of electron neutrinos is created in the 
interior of the sun, passing through 10 11 gm cm -2 of matter and eventually reaching a terres- 
trial detector located at a distance of 10 s km from the sun. The tests discussed in this paper 
are independent of solar models and are made possible by the fact that low energy (MeV) 
nuclear fusion reactions produce only electron type neutrinos. For neutrino squared mass 
differences less than 10~ 4 eV 2 , the solar neutrino tests are more sensitive than laboratory 
tests P| of lepton flavor conservation. 

We consider measurements of: 1) the energy spectrum of recoil electrons in charged 
current (absorption) reactions 2) the ratio of the number of charged-to- neutral current 
events 0; and 3) the combined measurement of the charged current energy spectrum and 
the charged to neutral current ratio. The shape measurement is sensitive to an energy- 
dependent depletion of the created flux of electron flavor neutrinos, and the neutral current 
to charged current comparison is sensitive to a non-zero conversion probability to a different 
(active) neutrino. 

How can we test lepton flavor conservation with solar neutrinos? The energy spectrum of 
8 B neutrinos is the same in the laboratory and in the sun, modulo negligible (gravitational 
redshift) corrections of C(10~ 5 ) ||. Fortunately, the spectrum, X(E V ), can be determined 
with relatively small uncertainties from laboratory data on the 8 B(/3 + ) 8 Be(2a) decay chain 
0. The measurement of the electron spectrum produced by neutrino absorption is therefore 
a test for new physics independent of complications related to solar physics. The ratio of 
neutral current to charged current events is also independent of uncertainties that affect the 
calculation of the total flux.[] All of the observed solar neutrinos must be of the electron 
type unless the separate conservation of electron flavor is violated. 

Neutrino oscillations are used in this paper to illustrate the potential effects of flavor 
transitions, but the considerations described here can be applied to other proposed electron- 
flavor non- conserving mechanisms, such as neutrino decay [0,U, non-standard electromag- 
netic properties P-P"T[|, neutrino violation of the equivalence principle [I2|, and supersym- 



metric flavor-changing neutral currents |T^JT^]. Many of the key papers and other relevant 
references are reprinted in [TE . 

Will the measurements with SNO be sensitive enough to prove — if it is present — that new 



x The calculated total flux of 8 B neutrinos at earth (all flavors) depends on solar physics and on 
the extrapolated low-energy cross section for the reaction 7 Be(p, 7) 8 B. 
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neutrino physics is occurring? Will the uncertainties (systematic and statistical) be suffi- 
ciently small to identify electron flavor non-conservation if it occurs as previously suggested? 
The answer is: "Yes, if SNO performs as expected" 

The SNO collaboration is completing the construction of a 1000 ton deuterium detector 
in the Creighton Mine (Walden, Canada) ||16| . The detector will measure the rates of the 
charged (CC) and neutral (NC) current reactions induced by solar neutrinos in deuterium: 

v e + d— > p + p + e~ (CC absorption) , (1) 

v x + d — > p + n + i> x (NC dissociation) , (2) 

including the determination of the electron recoil energy in Eq. ([]]). Only the more energetic 
8 B solar neutrinos will be detected^ since the expected SNO threshold for CC events is an 
electron kinetic energy of about 5 MeV and the physical threshold for NC dissociation is the 
binding energy of the deuteron, = 2.225 MeV. 

Neglecting all systematic uncertainties, some previous authors |],[TB|-f2r]] have considered 
how well the tests of flavor conservation by SNO can discriminate between new neutrino 
physics scenarios and standard model expectations. The most explicit discussions are given 
in [[19] and [p0|| , which constitute especially good introductions to the subject. We evaluate 
the effects of systematic uncertainties, theoretical and experimental, on the discriminatory 
power of the SNO tests for new physics. We consider uncertainties related to the laboratory 
shape of the neutrino energy spectrum, the calculated cross sections for charged-current and 
neutral-current reactions with deuterium, the energy calibration and resolution, detection 
efficiencies, and the CC detection threshold of the SNO detector. 

The primary goal of this paper is to refine the best-estimates and uncertainties for the 
theoretical ingredients that will determine how powerfully SNO will test electron flavor 
conservation. In addition, we carry out a preliminary overall estimate of the sensitivity of 
the detector to different types of oscillation phenomena, including realistic estimates of the 
experimental characteristics and their uncertainties. The various backgrounds |l] will not 
be known until after the detector is operating and are not included here. 

Our analysis is not a substitute for the detailed Monte Carlo simulations of the operating 
detector which will be performed by the SNO team. The discriminatory power of the detector 
will be established definitively by the simulations to be performed by the SNO collaboration, 
which will include all the theoretical ingredients discussed here, the detector elements that 
we high-light, and other aspects of the detector (such as the backgrounds) that will be 
determined during the operation of the experiment. Our calculations can, however, be a 
useful guide as to what is likely to be possible and what uncertainties are most important to 
try to reduce. We note that the SNO collaboration has been working for a number of years to 
develop the detector and calibration techniques in ways that will minimize the experimental 
uncertainties. 

We shall show that the recognized systematic uncertainties permit the observation of 
new physics at SNO, but the systematic errors may well dominate the total uncertainties 



The contribution of hep neutrinos [17| is negligible and will be discussed in Sec. II A. 
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after a relatively short exposure (~ 1 year) to solar neutrinos. Our analysis can be extended 
easily to include additional sources of uncertainties. 

We concentrate here on the most direct tests for new physics, which involve the shape of 
the neutrino spectrum and the charged-current to neutral-current ratio. If SNO does find 
evidence for new physics, the next step will be to discriminate among competing models of 
new physics. Important information will be provided by the time-dependence of the observed 
solar neutrino signal (day-night and seasonal variations) [[TJ. We do not address questions 
related to the time-dependence in this paper. 

The SNO collaboration plans an overall test of the detector by measuring the energy 
spectrum of an intense 8 Li beta-decay source to be placed in the SNO detector [^TJ. We 
include in Appendix A our determination, using the best-available data, of the 8 Li(/3) spec- 
trum, together with its estimated uncertainties. We also discuss some possible strategies for 
the 8 Li-test. 

This paper has the following structure. In Sec. [IT], we describe the neutrino-related 
ingredients of our calculation: the 8 B neutrino energy spectrum, and the charged and neutral 
current neutrino cross sections for deuterium. In Sec. [IT]] , we discuss the detector-related 
ingredients: the energy resolution, the absolute energy scale, the detection efficiencies, and 
the CC energy threshold. We use in Sec. [TV] the neutrino-related and the detector-related 
ingredients to calculate the flavor-conserving expectations for the shape of the CC electron 
recoil energy spectrum and the CC/NC event ratio; we include realistic estimates of the 
likely uncertainties and the correlations among the uncertainties. In Sec. [V], we calculate 
the effects upon measurable quantities of representative neutrino oscillation scenarios, and 
assess quantitatively the statistical significance with which new physics might be observed. 
We summarize our work in Sec. [VI[ Appendix A presents a calculation and discussion 
of the 8 Li(/?) spectrum and its use as a test of the overall performance of the detector. 
Appendix B discusses the extent to which the average value of the electron recoil energy is 
a good estimator of possible deviations from the CC shape expected in the absence of flavor 
violations. 



II. NEUTRINO-RELATED INGREDIENTS 

In this Section, we discuss the neutrino-related ingredients of the analysis that have 
appreciable, recognized uncertainties. These ingredients are: the 8 B neutrino spectrum 
(Sec. |11 A[ ), the charged-current absorption cross section (Sec. |11 B| ), and the neutral-current 



dissociation cross section (Sec. |II C| ) . We discuss the detector-related ingredients in the 
following section. 



A. 8 B neutrino spectrum 



The only component of the solar neutrino spectrum |17| that is important for the SNO 
experiment is the 8 B spectrum, X(E U ). A derivation of the best-estimate 8 B spectrum, X(E U ), 
along with the maximum allowed deviations, X ± (E U ) (±3 effective standard deviations away 
from the best-estimated spectrum) is presented in reference H . 
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The 8 B neutrinos are produced in the decay 8 B(/3) 8 Be followed by 8 Be(2a) disintegration. 
The broad 8 Be intermediate state is responsible for important deviations of the neutrino 
spectrum A from the usual allowed shape. The population of the 8 Be state is determined 
experimentally by measuring the delayed a-decay spectrum. The absolute energy calibration 
of the measured alpha spectrum is the main systematic error. 

For the calculation of the 8 B neutrino spectrum, the experimental and theoretical un- 
certainties can be included in a single effective alpha-energy offset, b: E a — > E a + b. 
The independently-measured positron spectrum in 8 B(/9) decay [fj2| provides a fundamental 
additional constraint that was used in || to choose a "best" reference a-spectrum and 
to bound its offset: b = 0.025 ± 0.104 MeV (±3cr uncertainties). An "infinitely-precise" 
measurement of the 8 B(/3) positron spectrum could reduce the effective ±3<r range of b to 
±0.075 MeV, where 0.075 MeV is the residual theoretical uncertainty. The uncertainties of 
the 8 B neutrino spectrum would be reduced in the same ratio. Since the uncertainties in 
the neutrino spectrum are a significant source of errors for the CC-shape test with SNO, a 
reduction of the allowed range of b through more precise measurements of the 8 B positron 
spectrum would be useful. 

The hep neutrinos |17| have a maximum energy of 18.8 MeV and could also in principle 
contribute to the neutrino spectrum observed by SNO. The calculated total flux of hep 
neutrinos is uncertain by a factor of about six p3] because of theoretical difficulties in 



estimating the low-energy production cross section. Using the nominal value <phep = 1 x 
10 +3 cm _2 s _1 given in f25| , we estimate that hep neutrinos contribute less than 0.07% 



of either the total NC or CC event rates. Therefore, hep neutrinos can be neglected in 
calculating the CC/NC ratio. Moreover, we have verified that the hep contribution to the 
high-energy tail of the spectrum is much smaller than the shape uncertainties estimated 
below in Sec. [IV Af 



B. Charged current ud cross section 

The cross section for the charged-current reaction (1) has been calculated a number 
of times in the last 30 years, since the original proposal by Jenkins [^6j to use charged- 
current capture on deuterium to measure the 8 B solar neutrino flux. Kubodera and Nozawa 
p7| have recently presented an insightful and thorough summary of the calculations of 
both the charged current and the neutral current cross sections. In this subsection, we 
assess the reliability of the theoretical calculations of the total and the differential CC cross 
section. We establish the robustness of the calculated cross sections, which is exemplified 
by the excellent agreement between the simple effective range calculations and the more 
sophisticated treatments. In addition, we stress the importance for determining the shape 
of the electron spectrum of including the final state interactions among the protons . We 
discuss the neutral current cross section in the following subsection. 



The kinematics of reaction (1) leads to the following expression |28] for the neutrino 
energy, E u : 

E u = Q + T e + — + {Pe ~ P ^ , (3) 

trip 4m p 
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where p v is the neutrino 3-momentum, T e and p e are the electron kinetic energy and 3- 
momentum, P is the relative momentum of the protons in the proton center-of-mass (c.m.s.) 
system, and the threshold Q = 1.442 MeV. The third term describes the kinetic energies of 
the two protons in the proton c.m.s. and is an important contribution to the energy balance. 
The fourth (last) term in Eq. (|3]) describes the small recoil energy of the two-proton center- 
of-mass system. 

The charged-current absorption, reaction (1), is described by the well-known electroweak 
Hamiltonian and by less well-known nuclear physics effects that can be treated at various 
levels of approximation. A neutrino energy of 10 MeV corresponds in natural units to 
(20 fm) _1 . Therefore, we expect that the details of deuteron nuclear physics will play only a 
minor role. This expectation is confirmed by comparing the s-wave calculations performed 
in the 1960's by Kelly and Uberall |25 and by Ellis and Bahcall |30[ using Bethe's effective 
range approximation |3l| with the recent sophisticated calculations by Ying, Haxton, and 
Henley [^] and by Kubodera and collaborators JE^TJ. The recent calculations include 
higher partial waves, relativistic effects, forbidden matrix elements, and exchange-currents. 

Figure |I]a shows the excellent agreement between the effective-range calculations and the 
more sophisticated treatments. In the figure, three independent estimates of the total CC 
cross section are compared: Ellis-Bahcall (EB) ||30|| , Kubodera- Nozawa (KN) |27| , and Ying- 
Haxton-Henley (YHH) |32j. The main difference between the calculated cross sections of KN 
and YHH is an energy-independent normalization factor (about ~ 6% uncertainty, as also 



estimated in ||27|| ). The EB normalization shows a slight energy dependence, that amounts to 
a ~ 4% additional variation over the important energy range of 5 MeV to 10 MeV. Figure |I]b 
(referring to the neutral current cross- section) will be discussed in Sec. [11 C| . 

When quoting the Ellis-Bahcall cross-section [3C], we use a slightly- improved calculation 
of the differential CC cross section in which we include the previously neglected p-p c.m.s. 
recoil term [the fourth term in Eq. (§)]. We have also used the more recent choices of 
parameters of the effective range approximation given in reference []34| , in order to obtain 
our best estimate of the differential cross section dace (E v , T e , cos 9) / dT e d cos 9 for any given 
value of E u , T e and of the electron scattering angle 9. 

Figure |2| shows the results of the improved Ellis-Bahcall calculation of the normalized 
differential cross section a^ c dace / dT e d cos 9 (dotted line) as a function of the dimensionless 
variable T e j[E v — Q), for representative values of E v and 9. The Kubodera- Nozawa results 
appear as a solid line in the same figure.^ 

The close agreement between the EB and KN normalized recoil electron spectra is strik- 
ing. The angular dependence calculated using the effective range and the Hamiltonian 
approximations is essentially identical. On the basis of Fig. ^ we conclude that the un- 
certainties associated with the (T e , 6 l )-shape of the normalized differential cross-section for 
the reaction (1) are much smaller than other recognized uncertainties. In practice, we 
parametrize the reference cross sections EB, KN, and YHH in the form: 



3 The extensive numerical tables of the differential charged current cross section calculated by 
Kubodera and Nozawa in [27| are not published. We thank the SNO collaboration for providing 
us with a computer-readable copy of these tables. 
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dacc(E v , T P 



dT P d cos 9 



X 



1 da cc (E UJ T e , cos 6) 
acc(Ev) dT e dcos9 



(4) 

EB 



for X = KN, YHH. The differences between the EB, KN and YHH cross sections are em- 
bedded in a multiplicative factor that depends exclusively on E v and not on the angular 
distribution. 

Contrary to what is sometimes assumed in the literature (see, e.g., p9|j35[1 ), the electron 



spectra in Fig. 0, although peaked, cannot be approximated well by delta functions in the 
electron energy. In other words, there is not a one-to-one relation between the incoming 
neutrino's energy and the energy of the electron that is produced. The final state in the 
charged current reaction cannot be approximated by a pure two-body state. Even as early as 
the seminal Kelly-Uberall calculation |29|, it was noted that the attractive 1 S p-p interaction 
is not sufficient to bind the protons as an effective single particle, because of the presence of 
the repulsive Coulomb force. The two-body approximation would be equivalent to omitting 
the third and fourth terms in Eq. (^). Removing only the fourth (the smallest) recoil term 
in Eq. (|3|) would cause the dotted cross sections in Fig. |2] to be systematically peaked at 
slightly higher electron energies (about +2% at E v = 12 MeV). The excellent agreement 
between our improved Ellis-Bahcall calculation and the Kubodera-Nozawa differential cross 
sections would be spoiled by omitting even this smallest term. 

C. Neutral current vd cross section 

We use the recent calculations of the neutral current cross sections (averaged over final 
states), <jnc(E u ), by Kubodera and Nozawa (KN) f2"7jl , and by Ying, Haxton and Henley 
(YHH) p2| . Only the total rate for reaction (fj), not the differential production rate as 
a function of energy, will be measured by SNO. However, the energy dependence of the 
cross section for the neutral current reaction is relevant for the SNO experiment, since the 
calculated differential cross section will be used in the SNO Monte Carlo simulations to 
model the production, and subsequent detection, of the neutrons produced by the neutral 
current reaction. 

Figure [l]b compares the KN and YHH neutral current cross section as a function of 
neutrino energy. The difference is ~ 6%, about the same magnitude and in the same direction 
as for the charged current cross section. There is a small residual energy-dependence below 
E v ~ 6 MeV. The difference shown in Figure [TJt> between the theoretical calculations is 
consistent with the theoretical error of ±10% that was estimated by Bahcall, Kubodera, 



and Nozawa |36| from various contributions (the impulse approximation, the nucleon-nucleon 
potential, meson exchange currents, and higher partial waves). 

The best estimate (and la uncertainty) for the neutral current cross section averaged 
over the 8 B neutrino spectrum is 

(a( 8 B)) = 0.478(1 ± 0.06) x 10~ 42 cm 2 . (5) 

We have chosen in Eq. ([5]) the la error of 6% to reflect the difference between the KN 
and YHH calculations. The standard solar model |25| prediction for the neutral current 
event rate due to 8 B neutrinos is 
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(</>a> = 3.2j™ SNU, (6) 

where the quoted la error in Eq. @ combines quadratically the uncertainties in the solar 
model prediction, the 8 B neutrino spectrum, and the neutral current cross section. The 
uncertainties in the solar model calculation dominate the error estimate. The event rate 
SNU is defined pTj as 10 -36 interactions per target atom (deuterium atom) per second. 



In the calculation of the CC/NC rate, the (already- small) theoretical cross section errors 
largely cancel and therefore do not affect significantly the ratio. As a default choice, we use 
the KN neutral current cross section in our calculations. The YHH cross section is used for 
comparison and to evaluate the theoretical uncertainties. 



III. DETECTOR-RELATED INGREDIENTS 

In this section, we discuss the detector-related ingredients of the analysis: the energy 
resolution, the absolute energy scale, the detector efficiencies, and the energy threshold for 
detecting CC events. Accurate determinations of all of these experimental quantities, and 
their associated uncertainties, are significant for the success of the SNO experiment. 



A. Energy Resolution 

The measured electron kinetic energy, T e , determined by SNO with the Cherenkov tech- 
nique, will be distributed around the true energy T' e with a width established by the photon 
statistics. 

The resolution function R(T' e , T e ) is expected to be well approximated by a normalized 
Gaussian, 



<r(Tl)y/2n 



exp 



2a(T ( 



(?) 



with an energy-dependent one-sigma width cr(T e ') given by |T7]. 



*{T' e ) 



0io1 



V 

e 



10 MeV 



(8) 



where 010 is the resolution width at T' e = 10 MeV. A plausible estimate of the parameter 
oio is 1.1 MeV (cxio = 1.8 MeV for Kamiokande, see f38fl), and a w itself may be uncertain 



by 10% 1 39|| . We will use in what follows cxio = 1.1 ± 0.11 MeV (la errors) as an illustrative 
estimate. 

In Figure §, solid line, we anticipate the results of our best-estimate of the standard 
shape of the electron energy spectrum (see Sec. |TV| ). The dotted line in Fig. |3| represents the 
same spectrum without the inclusion of the energy resolution. The area under the curve is 
normalized to unity in both cases. 
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B. The Absolute Energy Scale 



What is the absolute accuracy of the SNO energy scale? How precisely will the average 
SNO energy measurement correspond to the true electron energy? 

The energy resolution function of the previous section describes how the measured kinetic 
energy, T e , is distributed around the true energy, T' e , assuming that the centroid of the 
distribution, T e ave , coincides with T' e . The calibration of the energy scale will be performed 
with a series of 7-ray sources, the most important of which are monoenergetic. The primary 
energy calibration source will be the 6.130 MeV 7-ray from the decay of the first 3~ excited 
state in 16 0. We define the systematic error in the absolute energy calibration, 5, by the 
relation 

8 = T e ^ e -T' e . (9) 



A reasonable la estimate 391 of 5 is: 



* = ±lookeV (lolfev) • °^ a£1 • (10) 

which corresponds to a ±1% error at 10 MeV. For comparison, the Kamiokande collaboration 
achieved |38|] a ±3% energy scale error. 



The case a = (a = 1) would correspond to an energy- independent scale shift (scale 
factor). The intermediate case a = \ would apply to a scale uncertainty dominated by the 

width of calibration lines (error oc Jt£)- 

In general, the phenomenological parameter a will depend both on the physical sources 
of the scale uncertainties and on the calibration technique. It may even remain an unknown 
parameter after calibration. However, as we will see, the energy-scale induced uncertainties 
of the SNO observables depend only mildly on a. The worst case appears to be a = (i.e., 
a uniform energy bias), which we adopt for a conservative estimate of the errors. 

In practice, we introduce the energy scale shift 5 by modifying the energy resolution 
function [Eq. (^)] with the replacement: 

R(T^T e )-^R(T^ + 5,T e ) (11) 



and 5 given by Eq. fllPp. The reader can verify that the transformation given in Eq. ( |iT|) 
is an appropriate representation of the energy-calibration uncertainty by writing the total 
rate for the process under consideration as a triple integral over the neutrino energy (E u ), 
the true electron recoil energy (T e '), and the measured electron recoil energy (T e , between a 
specified minimum and maximum value). 



C. Detection efficiencies 

The detection efficiencies, ecc an d 6nc> f° r detecting charged and neutral current events, 
will be measured with calibration experiments at SNO along with their uncertainties, cr(ecc) 
and cr(eNc)- The calibration experiments will also measure possible variations of the CC 
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efficiency with the (true) electron recoil energy, ecc = e cc(T^), which can in principle affect 
the CC-shape measurements. If the experiments work as expected f40" l, then cx(ecc) < o"(cnc) 
and 



a(e NC ) ^ 2% (la) . (12) 

In Sec. |Vj, we shall include the efficiencies ecc^e) and e^c in the general expressions 
for the predicted quantities. As default values, we will assume that ecc(T e ) is constant and 
equal to 1, and that 6nc = 0.50±0.01 (la error). It will be shown in Sec. [IV C| that plausible 
energy variations of ecc induce deviations in the CC spectrum that are much smaller than 
other sources of error. 

In any event, after ecc(T' e ) and £nc are measured in the SNO detector, their effects on 
the predictions, and their uncertainties, can be included easily using the formalism given 
here. 



D. Threshold Energy 



The measured kinetic energy threshold, T min , for counting charged current events is 
expected to be fixed around 5 MeV [l|]. Below ~ 5 MeV, the signal-to-background ratio is 
expected to decrease very rapidly. In principle, one would like to have the threshold as low 
as possible in order to increase the number of events that are detected for a given exposure 
and in order to observe more of the curvature of the spectrum at lower energies (cf. Fig. 
The actual background level in the operating SNO detector will determine how low the 
energy threshold may be set. In Sec. |V CJ , we will show that the discriminatory power of the 
experiment is not very sensitive to the actual threshold level as long as the threshold is in 
the vicinity (±1 MeV) of the nominal value, T min = 5 MeV. 



IV. STANDARD NEUTRINO PHYSICS AT SNO 

In this section, we calculate the standard predictions, and their associated uncertainties, 
for the shape of the CC electron recoil energy spectrum and the ratio of total number of CC 
to NC events. We assume standard neutrino properties (lepton flavor conservation and zero 
neutrino masses) and the ingredients discussed in the previous sections (Sees. |Tj and p|) . We 



adopt the Kubodera-Nozawa (KN) neutrino interaction cross sections as "standard," since 
they are the most recent and complete for both the charged and the neutral current reac- 
tions. In particular, we evaluate the standard CC differential cross-section as indicated by 
Eq. (|4]) with X = KN, i.e., with absolute normalizations given by the Kubodera-Nozawa p7[ 



calculations and relative differential cross sections given by the Ellis-Bahcall |3(J calculation 



The shape of the recoil electron energy spectrum is given by the normalized distribution 
of charged current events (Vcc) as a function of the measured kinetic energy, T e : 

^da C c 



[ dE u \(E u ) f dT^R(T>, T e )e cc (T>) 

dTL , (13) 



1 dN cc 

Ncc dT e f Jm f ,„ WC1 , f ^dacc 



dT e J dE v \{E v ) J dT' e ^R{T' e , T e )e cc (Z 
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where T min is threshold for the measured electron kinetic energy and eccC^) is the efficiency 
for detecting an electron of true energy T' e . The energy resolution function, R{T' el T e ), is given 
by Eq. ([?]), with an allowance for an energy scale shift [Eqs. (|T0| , [XT])] . The CC differential 



cross section, %pr, is implicitly integrated over the entire solid angle since events will be 
detected for all efectron recoil angles. 

The ratio of charged to neutral current events may be written: 



Ni 



cc 



/ dT e [ dE u X(E u ) f dT' e ^R(Tl T e )e cc (Z) 



Nnc e NC J dE v \{E v )a nc {E v ) 



where £nc is the overall efficiency of neutral current event detection. 

It is necessary to adopt specific values for the efficiencies ecc and £nc in order to evaluate 
the relative number of CC and NC events and thus the statistical errors. We adopt plausible 
default values, ecc = 1 and €nc = 0.50. However, after calculating the statistical and 
efficiency errors, we prefer to convert the results to, and to quote, an essentially efficiency- 
independent CC/NC ratio, _Rcc/-Rnc> defined as 

-RcC Cnc ^CC / 1c -x 

(15) 



-Rnc ecc N] 



NC 



In writing Eq. flT5p, we have assumed that ecc is equal to a constant, which turns out to be 
a good approximation. 



A. Standard model predictions for the CC spectrum shape 

What effects do the different uncertainties have on the shape of the electron energy 
spectrum? Figure [| answers this question by showing the standard spectrum (solid line) 
and the effective 3a shape errors (dashed lines) due to neutrino-related and to detector- 
related uncertainties. We postpone the discussion of the sensitivity of the results to the 
assumed CC threshold energy to Section [VQ (see especially Tables |T| and |TV|). 



The dominant neutrino-related uncertainties are due to the 8 B neutrino energy spectrum 
\{E V ) (Fig. ||]a). The theoretical CC cross section uncertainties are very small; Fig f|b shows 
the "greatest" deviation, induced by the use of the EB instead of the KN cross-sections. 

The detector related uncertainties are due to statistics (Fig. |c), energy resolution 
(Fig. |]d), and energy scale (Figs. |]e,f). The statistical errors bars in Fig. ||c refer to 
a hypothetical sample of 5000 CC events collected above threshold and divided in 10 
bins. In Fig. §d, the dotted curves have been obtained by using energy resolution widths 
(Tio = 1.1 ± 0.33 MeV (±3cr errors, see Sec. [Ill A|) . The last two subplots of Fig. | show 
the effect of the absolute scale uncertainty, with a = or a = 1 (see Sec. |IHB| ). The two 
cases are almost indistinguishable. We adopt in the following the value a = 0, since it gives 
slightly more conservative error estimates. 

Figure ^ shows that the systematic errors due to the adopted uncertainties in the energy 
scale, in the resolution width and in the 8 B neutrino spectrum are at least as important 
as the statistical errors, with the additional complication that the non-statistical errors are 
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correlated point-by-point. The correlations of the uncertainties imply that the analysis of 
a realistic spectrum divided in N bins could become rather cumbersome, requiring con- 
sideration of an N x iV covariance matrix with large off-diagonal elements for the usual x 2 
statistics. Moreover, the x 2 test is not powerful when bins are affected by significant system- 
atic errors, since the additional information embedded in sequences of equal-sign deviations 
(typically all positive or all negative in one half of the spectrum, as in Fig. [|) is lost. 

The simplest quantity that characterizes a generic electron spectrum, while avoiding the 
use of bins, is the average value of the measured recoil energy, (T e ) (see also |T9| , pO| ) .f\ The 
basic question then becomes: "Can SNO detect significant deviations of (T e ) measured from 

(Te) standard • 

We evaluate (T e ) with the aid of Eq. (|13|) by using the definition 

(T e }=[ dT e T e — ^. (16) 
V ' JT min e N cc dT e 1 J 

We determine the effects of uncertainties in different ingredients by carrying out the 
integration indicated in Eq. ([IB]) with different assumptions. For example, we estimate the 
la uncertainty associated with the neutrino spectrum by evaluating (T e ) using the neutrino 
spectra X + (E U ), X~(E U ), which are ±3cx away from the best-estimate neutrino spectrum 
||. Then the la difference is a((T e )) = |[(T e (A + )) — (T e (A - ))]. Analogously, the la errors 
due to the energy resolution uncertainties were estimated by recalculating the spectra, and 
thus (T e ), with a io = 1.1 ± 0.33 MeV (±3<r), and dividing the total shift by six. A similar 
procedure was adopted for determining the energy scale error. For the CC cross section 
uncertainty, we have attached a la significance to the deviation obtained when the Ellis- 
Bahcall CC cross sections were used instead of the Kubodera-Nozawa CC cross sections, 
i.e., a((T e )) = (T e (EB)> - (T e (KN)>. 

Our standard estimate, for T min = 5 MeV, is then: 

(T e ) = 7.658 ± 0.025 a ± 0.011 6 ± 0.029 c ± 0.024 d ± 0.052 e MeV , 

= 7.658 • (1 ± 0.009) MeV , (17) 

where the errors (±lo~) are due to: a) statistics of 5000 CC events; b) difference between 
EB and KN cross sections; c) uncertainties in the neutrino spectrum; d) energy resolution; 
and e) the absolute energy calibration. The statistical error of the average is calculated, 
according to the central limit theorem, as <T sta t = yVar/Vcc, where Var is the variance of 
the standard distribution above threshold, Var = (1.74 MeV) 2 . 

The assumed uncertainty in the absolute energy scale, Eq. fllPf) , dominates the total 
error. 

In the best-estimate calculation, the detection efficiency ecc has been taken constant. 
A linear dependence of ecc on T e (neglecting temporarily the distinction between measured 
and true energies, T e and T' e ) would modify the CC spectrum as: 



4 For a discussion of the extent to which the average electron recoil energy is a good estimator of 
the deviations from the expected shape in the absence of electron flavor violation, see Appendix B 



and Ref. 141 
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1 dN cc 1 dNco( 1 + f ,^pp\ (lg) 



N cc dT e N cc dT e \ (T e ) 

and the average value as 

(Te) — (T e ) (l + (19) 

where (3 is the slope of ecc — e cc(T e ). A plausible variation (or uncertainty) of ecc over 
the interval 5-15 MeV is a few percent, say 3% for definiteness. This would correspond 
to /3 — 0.023 and to a 0.12% variation of (T e ), comparable to the smallest error shown in 
Eq. dT7|). We conclude that uncertainties in the CC efficiency will not be an important 
source of error, if SNO performs as expected. 

The characterization, Eq. (|i~6l), of the spectrum given in Eq. (fl3|) is the first (non-trivial) 
step in a complete description by a series of moments: the zeroth moment (the area), the 
first moment (the average), the second moment (the variance), and higher-order moments. 
In the present case, the zeroth moment is equal to unity by definition, and small varia- 
tions in the shape of the spectrum affect primarily the average value. More sophisticated 
unbinned tests, such as the Kolmogorov-Smirnov (K-S) test, may be useful to apply after 
the SNO collaboration has estimated by empirical calibrations the systematic errors in the 
experimental input quantities. One could then determine by Monte Carlo simulations the 
distribution function for the K-S statistic with an inferred model for the systematic errors. 



B. Standard model predictions for the CC/NC ratio 



In Eq. ([TBI), we have defined a charged to neutral current ratio, Rqc/Rnc, which is 
independent of the average absolute values of the efficiences but incorporates the efficiency- 
induced errors. 

The calculation of the standard value and ±ler errors for Rqc/Rnc is done with the 
same logic as for (T e ). The final result is: 



Rcc 



1.882 ± 0.058 a ± 0.010 6 ± 0.008 c ± 0.009 d ± 0.034 e ± 0.038 / 



-Rnc 

= 1.882(1 ±0.042) . (20) 

The individual contributions result from: a) statistics of 5000 CC events and 1354 NC events 
(cnc — 0.5); b) difference between YHH and KN cross sections; c) neutrino spectrum; d) 
energy resolution; e) energy scale; f) NC efficiency. 



C. Correlation of errors 

Some of the systematic errors affecting (T e ) and Rqc/Rnc have the same origin and are 
correlated. In particular, a variation of the neutrino spectrum X(E U ) causes both (T e ) and 
i?cc/i?NC to increase or decrease at the same time. A variation of cr 1Q , the energy resolution 
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width, produces instead opposite effects (anti-correlation). Energy scale errors are positively 
correlated. 

In Table |, we summarize the separate error components and their correlations. The 
correlation of the total errors, p = 0.32, is small because of accidental cancellations, but it is 
not entirely negligible. When the actual error budget for SNO is determined experimentally, 
the approximate cancellation of the correlations may not be as strong. 



V. COULD SNO PROVE THE OCCURRENCE OF NEW PHYSICS? 

In the previous section, we calculated the standard predictions for the CC-shape and the 
CC/NC ratio, (T e ) and -Rcc/-RnCj along with a realistic estimate of these uncertainties. In 
this section, we show that the anticipated uncertainties are sufficiently small to allow SNO 
to prove the occurrence of new physics with a high degree of confidence. 



A. Neutrino oscillations 

We examine the implications of three representative scenarios in which the solar neutrino 
problem is solved by neutrino oscillations: the two (best-fit) Mikheyev-Smirnov-Wolfenstein 
(MSW) solutions at small and at large mixing angle (SMA and LMA), and the purely 



vacuum (VAC) oscillation [43] solution (see [44] and refs. therein). 



Figure [5] illustrates some of the principal differences among the three oscillation scenarios 
and also shows that all of the oscillation solutions differ significantly from the standard 
model expectations (STD). The survival probabilities for electron- type neutrinos vary greatly 
among the oscillation scenarios, as can be seen clearly in Fig. |^a. The energy spectrum 
of electron-type neutrinos at earth, shown in Fig. |5|b, is affected strongly by oscillations. 
Figure |5]c represents the different CC electron recoil spectra that are predicted for SNO, 
with the normalization: Area = 1. Until the neutral current is measured, the different 
recoil spectra must be compared with the same normalization, as in Fig. |5|c, because we do 
not know a priori the total number of electron-type neutrinos that are created in the solar 
interior. Once the neutral current is measured, we can compare the different oscillation 
and no-oscillation scenarios in a more informative way. Figure |5]d makes use of assumed 
measurements of the total neutral and charged current rates and shows the standard electron 
spectra normalized to Area = iVcc/^NC-Q A comparison of Figs. |^c and Fig. |5]d makes clear 



5 The best-fit mass and mixing values (Am 2 , sin 2 26) are taken from [Q: (5.4 x 10 6 eV 2 , 7.9 x 
10~ 3 ) for small angle MSW, (1.7 x 10~ 5 eV 2 , 0.69) for large angle MSW, and (6.0 x 10~ n eV 2 , 
0.96) for vacuum oscillations. Numerical tables of the oscillation probabilities for these 
best-fit scenarios were prepared by P. Krastev 44] and are available at the following URL: 
http:// www . sns . ias . edu/~ j nbj . 



6 If we had normalized the area to the efficiency-independent ratio, Area = Rqc/ R^Ci then the 
scale of the ordinate in Fig. 5d would be decreased by a factor of 2, but the relative shapes would 
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the importance of the neutral current measurement for interpreting the shape of the electron 
recoil spectrum. 

Figure ^| shows the ratio of the normalized electron spectra after oscillation (as displayed 
in Fig. [|c) to the standard spectrum. These ratios of spectra are approximately linear 
in energy. Therefore, the main effect of oscillations on the normalized CC spectrum is to 



change the first moment of the energy distribution, i.e., the mean value (T e ). In Ref. pO 
it is shown that the approximate linearity of the ratios of recoil energy spectra is a general 
feature of the resonant MSW effect. The representation of the spectral information by one 
parameter, (T e ), is efficient because the ratios of spectral shapes are approximately linear 
(i.e., are determined rather well by just one parameter). 



B. Statistical analysis 

For the three neutrino oscillation scenarios considered in the previous section, we calcu- 
late the observables (T e ), Rqc/Rnc, and their distance — in units of standard deviations — 
from the standard predictions, Eqs. flTT| ) and (18). 

For each test, the distance is defined simply as: 

N{a) = (X - X s t a ndard)/Cx ,standard j 

(21) 

where X = (T e ), Rqc/Rnc- In the combined tests, we have calculated the \ 2 including the 
correlation of the total errors (see Table |), and defined AA(cr) = yfx 2 ||. When N(o) ^> 3, 
the physical interpretation is that the statistical probability of the result under consideration 
is negligibly small; the normal distribution presumably does not describe the extreme tails 
of the probability distribution. 

The expected deviations for different oscillation scenarios are shown in Table |I[ 

The measurement of the Rcc / -Rnc ratios is a powerful test for occurrence of new physics; 
the three oscillation cases are each separated from the standard expectations by a distance 
that is formally more than 15cr (cf. comment above regarding M{a) ^> 3 ). 

The CC-shape test is less powerful. This is of course expected for the large-angle MSW 
case, but is somewhat surprising for the small-angle MSW case, that was generally expected 
to be separated from the standard model expectations at a high confidence level |119p0| . 
The reason that the significance level found here for the measurement of the CC shape is 
much less than previously calculated is that we have included estimates for the systematic 
uncertainties. The systematic errors in measuring the CC-shape may be twice as large as 
typical statistical errors (5000 CC events in our case), as evidenced in Table |. The statistical 
power of the combined tests (CC-shape and NC/CC) is dominated by the measurement of the 
CC/NC ratio. The effect of the correlation of the errors is small but not entirely negligible. 

Figures [7] and |8] display graphically the information contained in Tables [I] and |T[ In 
Fig. ^ we show the standard predictions for (T e ) (upper panel) and Rqc/Rnc (lower panel), 
together with the separate and combined 3a errors. The values of (T e ) and Rqc/Rnc for 



remain the same. 
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the different oscillation channels are also displayed. In Fig. 7a, the efficiency error (labeled 
by a question mark) should be negligible if SNO works as expected. 

In Fig. |8| we show the results of the combined tests (correlations included) in terms of 
iso-sigma contours in the plane ((T e ), Rcc/Rnc), where M{a) = ^/x 2 . The three oscilla- 
tion scenarios can be well separated from the standard case, but the vertical separation 
(Rcc/Rtssc) is larger and dominating with respect to the horizontal separation ((T e )). 

The error bars on the SMA point in Fig. [7] and Fig. [S] represent the range of values 
allowed at 95% CL. by a fit of the oscillation predictions to the four operating solar neutrino 
experiments they are intended to indicate the effect of the likely range of the allowed 
oscillation parameters. 

The choice of (T e ) as a characterization of the CC-shape is not unique. (T e ) has been 
chosen because it is a single and well-defined number (the first moment of the electron 
distribution), whose systematic uncertainties can be determined independent of the event 
binning. If the measured electron distribution at SNO has significant deviations of the 
second moment (the variance), or higher moments, from the standard expectations, then 
(T e ) may not be the optimal statistical estimator. If a more sophisticated statistical test is 
used to test for non-standard curvature in the spectrum, then the assessment of statistical 
significance will require a full Monte Carlo simulation of SNO detector, with systematic 
effects calculated by brute force. 



C. Threshold Dependences 

All of the previous calculations were carried out assuming an recoil energy threshold, 
^min, of 5 MeV. The actual value of T min that will be used will depend upon the observed 
or estimated backgrounds in the operating SNO detector. 



Table |TJ summarizes the dependence of (T e ) and Rqc / Rcc on the adopted energy thresh- 
old. We give values for the standard model and for the three exemplary oscillation scenarios. 

In Table [TV] we give the distances of the oscillation scenarios from the standard predic- 
tions, in units of standard deviations. 

We see from Tables [TIT] and [IV] that the differences resulting from changing the threshold 
by ±1 MeV are not expected to be decisive for the SNO discovery potential. However, the 
diagnostic power of the measurement of the shape of the electron recoil energy spectrum 
would be significantly enhanced if the energy threshold were lowered. The small mixing 
angle MSW solution is 3.8a away from the standard model prediction if the threshold is 
4 MeV but is only 2.4a away if the threshold is 6 MeV. 



VI. SUMMARY AND CONCLUSIONS 

The Sudbury Neutrino Observatory has the potential to reveal new phenomena with a 
high level of confidence, but the detector must work well in order to discriminate among 
different physics options. The accurate calibration of the absolute energy scale, the energy 
resolution width for electron detection, and a high sensitivity for neutral current detection, 
are especially important. 
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We have determined both the the best-available estimates and the uncertainties 
of three important neutrino-related input quantities (available at the following URL: 
[littp:/ / www.sns.ias.edu/~jnpp that will be needed in the analysis of the SNO data: the 
laboratory shape of the 8 B neutrino energy spectrum, the charged current neutrino absorp- 
tion cross section, and the neutral current dissociation cross section. 

We have also estimated the effects on the tests of electron flavor violation of five detector- 
related aspects: the energy resolution, the absolute energy scale, the energy threshold, and 
the detection efficiencies of the charged current events and of the neutral current events. 

The principal uncertainties that affect the predictions are shown in Table [I] and in Fig [7]. 
For the measurement of the shape of the electron recoil energy spectrum, the largest esti- 
mated error is contributed by the uncertainty in the absolute energy scale, with significant 
additional errors arising from the energy resolution and from the shape of the 8 B neutrino 
energy spectrum. 

Unfortunately, the systematic uncertainties reduce the power of SNO to detect new 
physics via the measurement of the shape of the recoil electron spectrum. For example, a 
previous analysis [pH[] , which considered only statistical errors, indicated that a 3a distinction 
between the standard model prediction and the small angle MSW solution would be possible 
with only 1800 CC events observed with the SNO observatory. This same statistical-only 
analysis suggests that 5000 CC events would give more than an 8a distinction. With our 
adopted estimates of the systematic uncertainties and 5000 CC events, we find that, instead 
of 8a, the standard model prediction and the best-fit small angle MSW solution differ by 
only 3.1a, as judged by (T e ). Even with zero statistical error, the difference between the 
standard model prediction and the small angle MSW solution would be only 3.3<r. Indeed, 
the systematic uncertainties begin to dominate the statistical uncertainties for this case after 
(less than) a year of operation. The statistical significance can be improved by ~ la for the 
SMA solution if the variance of the spectrum as well as (T e ) is measured ||41]| . Fortunately, 
the shape of the recoil spectrum predicted by vacuum neutrino oscillations is distinctive and 
the difference between standard model physics and vacuum oscillations represents, with our 
adopted uncertainties, a 10cx distinction in the SNO detector. 

The main information content of the measured shape of the electron recoil spectrum 
can be summarized by evaluating the average recoil electron energy. For standard model 
physics, we find that the average electron kinetic energy is (T e ) = 7.658(1 ± 0.009) MeV, la 
total errors. 

Figure [5|c compares the normalized electron recoil spectra computed for the standard 
model, the LMA, the SMA, and the vacuum neutrino oscillation scenarios. The differences 
in the positions of the peaks of the spectra shown in Fig. ^c are larger than the differences 
in (T e ) given in Table |I| and Table |II|. Although previous authors have shown similar 
figures with error bars due only to the statistical fluctuations assigned to individual bins, 
we refrain from showing error bars in Fig. [|c since the systematic errors will likely dominate 
the statistical uncertainties and since systematic errors are correlated from bin to bin. 

There is not a one-to-one relation between the incoming neutrino energy in the charged 
current reaction and the energy of the electron that is produced. Figure 2 shows, in fact, 
that there is a significant spread in electron recoil energies for a specified neutrino energy. 
This result, unfortunately, contradicts the assumption of a one-to-one energy relation used 
by a number of authors [0,0 in describing potential applications of SNO measurements. 



17 



The neutral current to charged current event ratio is a sensitive probe of lepton fla- 
vor violation. For standard model physics, we find a charged-to-neutral current ratio 
Rcc/Rnc = 1.882(1 ± 0.042), la total error. 

The measurement of the absolute neutral current rate will test directly the solar model 
prediction |2§] of 3.2±g;| SNU, la total error, for the 8 B neutrino flux. This test of solar 
models is independent of uncertainties in the fundamental physics related to oscillations into 
active neutrinos. 

The predictions of the three favored oscillation solutions considered here (small mixing 
and large mixing MSW solutions, and vacuum oscillations) are all separated by more than 
16a from the predictions of the standard model with no lepton flavor violations, as shown in 
Table |T| and Fig. || The combined test, shape of the electron recoil spectrum and ratio of 
neutral current events to charged events, is, in our simulations, only slightly more powerful 
than the neutral current to charged current ratio alone. 

An accurate measurement of the neutral current rate is essential in order to exploit the 
full potential for new physics of the Sudbury Neutrino Observatory. 

Are the conclusions about the statistical significance of the flavor tests robust with respect 
to the SNO charged- current energy threshold? This question is answered in Table |T| and 
Table [TV]. These two tables show that the discriminatory ability of the SNO detector does 
not depend critically upon the detection threshold for the CC reaction. However, lowering 
the threshold to 4 MeV would separate by the shape measurement alone the small mixing 
angle MSW solution by 3.8a from the standard prediction, instead of the 3.1cr that applies 
for a threshold of 5 MeV. 

The relative insensitivity of the diagnostic power of SNO to the CC energy threshold 
suggests one possible strategy for dealing, especially in the initial stages of the experiment, 
with the most troubling backgrounds. Without seriously affecting the discriminatory power 
of the CC to NC ratio (see Table |T| and Table [TV]) , the CC energy threshold for events 
being analyzed can be set sufficiently high, at 6 MeV or perhaps even at 7 MeV, that no 
significant background contamination is plausible. 

Are the absolute event rates for the CC events sensitive to the assumed CC threshold? 
Table |V| gives the expected event rates for different assumed thresholds and neutrino oscilla- 
tion scenarios. The range of expected event rates is about a factor of two for CC thresholds 
from 4 MeV to 7 MeV for the standard model and for the SMA and LMA MSW solutions; 
the variation is about 40% for the vacuum oscillations. 

We have also varied the assumed value of the energy resolution width at 10 MeV, a w . 
We find that a one-third worsening of the energy resolution a w , which is defined by Eq. (|J), 
from the current best-guess of 1.1 MeV to 1.5 MeV decreases the difference between the 
standard model value for (T e ) and the small mixing angle (MSW) value from 3.1a to 2.8<r. 
Thus it is important that the energy resolution width is kept as low as possible. 

The expected background level in SNO decreases strongly with increasing energy [[[]. 
Increasing the electron energy threshold for the charged current reaction, can decrease the 
fractional contribution of background events. 

The principal lesson of Table |T| and Table |V| is that the SNO experiment will provide a 
powerful diagnostic for new physics even if the observed backgrounds are somewhat higher 
than expected. 

We have varied the mass and mixing parameters of the small mixing angle (MSW) case 
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within the 95% CL. limits of the fit to the four operating experiments p4 |. The difference 
between the standard value of (T e ) and the value calculated for the small mixing angle 
solution varies between 3.1 ±0.6 standard deviations, depending on which values one adopts 
within the allowed MSW region. (Somewhat more powerful discrimination can be achieved 
if both the dispersion and the mean recoil energy are calculated |^T|.) The formal difference 



between the standard and the MSW value of Rcc/Rnc is always within 15.7 ± 1.6 standard 
deviations at 95% CL. Thus the MSW solution at small mixing angles (SMA) can be tested 
with a high level of confidence with the SNO experiment, although the information coming 
from the CC-shape might not be sufficient by itself. 

We discuss in Appendix A the information that can be gained from an overall test of the 
SNO detector using an intense 8 Li(/5 _ ) source. We calculate a theoretical 8 Li(/5~) spectrum 
and compare with the available data. We conclude that the existing data for the 8 Li(/3~) 
spectrum are not sufficient to permit an accurate test of the SNO detector and that a new, 
laboratory experiment is required. 

Appendix B discusses the extent to which (T e ) is a good statistical estimator of possible 
deviations in the CC electron spectrum. 

Finally, we must ask: How general are the conclusions given in this paper? The method of 
analysis and the discussion of the principal ingredients and their uncertainties will be of use 
in considering how well SNO can test potential new physics scenarios. We have evaluated 
the sensitivity of the likely operation of the Sudbury Neutrino Observatory to the MSW 
and vacuum neutrino oscillations that best-fit (see J44|) the results of the four pioneering 



solar neutrino experiments. The correct physical explanation may differ from the oscillation 
solutions considered here. The "true" solution of the solar neutrino problems may involve, 
for example, a distortion of the CC electron recoil spectrum that is much more drastic than 
is implied by the oscillation solutions considered here. In this case, the shape of the electron 
recoil spectrum from the CC reaction might indicate new physics that is not apparent by 
comparing the rates of the CC and the NC reactions [|35|] . The analysis presented in this 



paper is illustrative of the power of the SNO detector, but specific, quantitative inferences 
depend upon what, if any, new physics exists in the accessible domain. 
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APPENDIX A: TESTING SNO BY MEASURING IN SITU 8 Li BETA-DECAY 

The SNO collaboration plans to perform an overall test of the experiment by measuring 
the /3-decay spectrum of an intense 8 Li source that will be placed in different locations 
in the detector. The measurement by SNO of the 8 Li(/3~) spectrum will be used as a 
demonstration that the results obtained for a known beta-decay spectrum are consistent 
with those measured in the laboratory. 

Figure [| compares our calculated (see below) 8 Li spectrum with the standard model 
electron spectrum from 8 B solar neutrino CC absorption [reaction (|])]. The test is based 
upon the fact that, despite the different physical processes that are involved in the two cases, 
the electron spectra from 8 Li /3-decay and from 8 B solar neutrino absorption on deuterium 
have somewhat similar shapes and cover essentially the same energy range (from to ~ 13 
MeV). The spectra displayed in Fig. ^|are separately normalized to unity above the standard 
SNO threshold of 5 MeV and do not include broadening due to the finite energy resolution 
in the detector (or other signatures of the SNO detector). 

We discuss in this appendix some of the things that can be learned from the 8 Li test. 
For background information, we first summarize in subsection [A 1| the relations between the 



8 Li and the 8 B electron spectra. We then describe in subsection |A 2j how a future precision 
laboratory measurement of the 8 Li spectrum could be used, in conjunction with a SNO 
measurement of the 8 Li spectrum, to help determine characteristics of the SNO detector. 

1. Relations Between the 8 Li and 8 B Spectra 

The 8 Li electron spectrum is produced by the beta decay 

8 Li -> 8 Be + e~ + v e . (Al) 

The CC electron spectrum, whose measurement is one of the primary goals of SNO, is 
produced by a two-step reaction (beta decay followed by neutrino capture): 

8 B -> 8 Be + e + + v e , (A2a) 

v e + d^p + p + e- . (A2b) 

All three of these reactions are different, but the intermediate 8 Be states are the same 
in both 8 B (reaction |A2a| ) and 8 Li (reaction |A1| ) decay. The 8 Be excited states are unstable 
and break up into two alpha particles. As a consequence, the shape of the 8 Li f}~ spectrum 
[Eq. ( |A2a )1 and of the 8 B j3 + spectrum [Eq. ([A2 bjj ] deviate significantly from the standard 



allowed shape. The shape of the B neutrino spectrum and its uncertainties are affected as 
well (see M and references therein). 

Measurements of the delayed a spectrum allow one to determine the profile of the inter- 
mediate 8 Be state and thus to calculate the deviations of the 8 Li electron, 8 B positron, and 
8 B neutrino spectra from their allowed shapes. The best-estimated 8 B positron and neutrino 
spectra and their uncertainties were discussed extensively in ||. 
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In Fig. p!0| a we show our best-estimate for the 8 Li(/? ) electron spectrum, together with 



its ±3<T uncertainties. The calculation that leads to Fig. |iO| a is similar to the calculation 
of the 8 B(/9 + ) spectrum performed in I|, modulo different values for the radiative correc- 
tions, forbidden corrections, and Coulomb effects. The form factors entering the forbidden 
corrections are, by isospin symmetry, the same as used for 8 B(/3 + ) ||. 

As discussed in Sec. II A, the 3cr uncertainties shown in Fig. |lO| a are related to a possible 
offset, b, in the energy of the a particles from 8 Be break-up: E a — > E a + b. The value of this 
offset affects, in particular, the calculated peaks of the @~, and v spectra in Eqs. (Al) 
and (A2). The effective 3a uncertainty of the offset b is estimated to be ±0.104 MeV ||; 
this estimate includes theoretical errors. 

There are at least three complementary experiments that could help to reduce the offset 
uncertainty, b, and thereby make the prediction of the 8 B solar neutrino spectrum more 
precise. The potential experiments are: 1) a high-precision measurement of the a spectrum 
from 8 Be break-up; 2) a high-precision measurement of the (5 + spectrum from 8 B decay; 
and 3) a high-precision measurement of the (3~ spectrum from 8 Li decay. In all three cases, 
dedicated laboratory experiments with a carefully calibrated spectrograph would be needed. 

The SNO detector is expected to have an uncertainty in the absolute energy calibration 
of ~ 100 keV (ler), and therefore probably cannot be used as a 8 Li spectrometer at the level 
of precision needed to further constrain b. 

2. The SNO Response Function: SsnoC^, T e ) 

In addition to an overall demonstration that the detector is working as expected, can one 
learn more about the characteristics of SNO by studying the Li(/3~) spectrum? The answer 
is: "Yes, provided that a precision measurement of the 8 Li(/?~) spectrum is made with a 
laboratory spectrograph." If future laboratory experiments determine accurately the 8 Li 
beta decay spectrum, then the measurement with SNO of this same spectrum can be used 
to constrain possible systematic effects that apply in the energy range that is also relevant 
for 8 B neutrino absorption []. 

Specifically, one could use the following strategy. Let A Ll (T e ') be the true lithium spec- 
trum as a function of the true electron energy, T' e . Suppose, as a first approximation, that 
this spectrum is known with "infinite" precision as a result of an error-free laboratory spec- 
trograph^ measurement: 

~ A Li . (A3) 

Then the lithium spectrum measured in the SNO detector, Ag^o; w ih be given by a convo- 
lution of A Ll (T e ') with the SNO response function, S^o{T' e) T e ), where T e is the measured 
electron kinetic energy and T' is the true electron energy. Thus 



7 The theoretical 8 Li spectrum in Fig. [To| a is affected by significant shape uncertainties (dotted 
lines) and is not a valid substitute for a high-precision laboratory measurement. 
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A 



Li 

SNO 



(T e ) = / A Li (T e ') ■ S SNO (T e ', T e ) dT> . (A4) 



In our notation [see Eq. (13)], the SNO response function can be written as the product 
of the energy resolution and the CC detection efficiency: 

SsnoCC T e ) = R(T' e , T e ) ■ e C c(T e ') . (A5) 

Therefore SsisroC^e, T e ) depends on three parameters: the energy resolution width a w 
[Eq. (8)], the absolute energy scale error 6 [Eq. (9)], and the slope of the CC efficiency 
function (3 [Eq. (18)]. (Of course, other free parameters could be eventually introduced to 
model Ssno more accurately.) Once A Ll and Ag^o are experimentally determined, one can use 
Eq. (A4) to fit the parameters of the SNO response function, Ssnch which describes the spec- 
tral distorsion effects induced by the detector. Since the SNO detector will also be calibrated 
with more traditional techniques (e.g., with 7 rays of known energy), the determination of 
the free parameters in Ssno will be overconstrained. The experimental overdetermination 
of Ssno will limit the effects of unknown systematic errors. Moreover, the comparison of 
the fitted parameter values with those estimated by Monte Carlo simulations will provide 
further consistency checks. 

In practice, one has to take account of uncertainties in the spectrum that is measured in 
the laboratory, Aj^ b , in order to infer the true spectrum A Ll . The corrections will depend on 
the response function of the laboratory spectrograph, S\^: 

AJtCTe) = / A Li (T e ') • S lab (T e ', T e ) dT> . (A6) 

The response function Si a b must be known with a precision higher than what one hopes 
to achieve for the overall SNO response function, Ssno- Any uncertainty, 5Sj a b, in the 
laboratory response function, will be propagated to Ssno- 

Are the available laboratory data on the Li(/3~) decay sufficiently good that their un- 
certainties would not introduce large errors in the SNO response function if determined via 
Eq. (|A~4T)? Unfortunately, the answer is "No." 



In Fig. |T0t b we show the published data from |45[] and from [jig] , which are superimposed 
on the theoretical spectrum. (We have taken the data at face value and made no attempt 
to deconvolve resolution effects.) The agreement of the data with themselves and with the 
theoretical spectrum is unsatisfactory. In Fig. |10|c, the j3 kinetic energies of each data set are 
linearly transformed, Tp —>■ ATp + B, so that the two (renormalized) experimental spectra 
match each other and fit the theoretical spectrum. A good fit-by-eye (Fig. |TU|c) is obtained 
with parameters: A = 0.94, B = 0.45 MeV for the data of reference [[45| and A = 1.06, 
B = —0.2 MeV for the data of reference [I5|. The results of these transformations suggest 
that true spectrum, A Ll , cannot be inferred from the available data with a precision better 
than a few hundred keV in the energy scale. 

A new, precision laboratory measurement of the 8 Li(/?~) spectrum is needed. 
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APPENDIX B: STATISTICS OF LINEAR DEVIATIONS FROM THE 

STANDARD SPECTRUM 



We state in the text that the average value of the electron kinetic energy, (T) , is a good 
statistical estimator of linear deformations of the recoil spectrum. What we mean by this 
claim is that (T) contains most of the information about spectral deformations, especially 
if the deformations are small. 

In this Appendix we show heuristically the basic correctness of the above statement in two 
simple but representative cases, purely statistical errors and a single dominant systematic 
error. More precisely, we show in Eq. (|B3|) (stat. error) and Eq. ( |E5| ) (syst. error) that, for a 
linear spectral deformation, the y 2 associated with deviations of (T) is approximately equal 
to the x 2 obtained by binning the observed spectrum in a histogram. However, if there are 
several comparable systematic errors, or if the deviation is nonlinear, significant additional 
information may be obtained from the higher moments [5T]. 



The arguments given below are in the spirit of a "physicist's proof" rather than a math- 
ematical theorem. We note that for probability distributions with long tails the analysis 
in terms of moments may not be appropriate. Fortunately, the spectrum of electron recoil 
energies does not have pathologically long tails so this last remark does not apply in the 
case we are considering. 

Let p(T) be the expected normalized electron recoil spectrum [JdT p(T) = 1], with av- 
erage kinetic energy (T) and variance a 2 . Let p'{T) be the observed normalized spectrum, 
with average energy (T)'. In the hypothesis of a perfectly linear spectral deformation, one 
can always write: 

P'(T) T- (T) 

where j3 is a slope parameter. Then the shift in the average energy is given by 

&(T) = (T)'-(T) = p-^, (B2) 
and its \ 2 statistic simply reads 

Xt^mt <B3) 



x G (T) 

where <J(t) is the total error affecting (T). 



8 Notice from Eq. (|B2| ) that A(T) and /3 are in one-to-one correspondence, so that a determination 
of the shift in the average energy (T) is equivalent to a determination of the slope of p'(T)/p(T) 
with the same fractional accuracy, and vice versa. The observables A(T) and (3 are interchangeable 
for a linear spectral distortion. We prefer A(T) because the average kinetic energy is well-defined 
also in the case of a non-linear distortion, while (3 is not. 
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For purely statistical errors, a 2 T ^ = cr 2 /N by the central limit theorem, where N is the 



total number of observed electrons. From Eqs. ( |B2|) and (B3) one has: 



cr 2 

x 2 = p 2 -^N (stat. dominated) . (B4) 
(T) 2 

A single small, purely systematic error (such as the uncertainty in the 8 B neutrino 
spectrum shape or in the absolute energy calibration) also produces, in first order, a linear 
deformation of the expected recoil spectrum. The spectral distortion can thus be represented 
as an uncertainty ap of the slope parameter (3. Then, from Eq. ( |B2"|) , the propagated error 
on (T) is (Ji T ) = crpcr 2 / (T), and one has from Eqs. ( |B2| ) and (|B~3|): 

xfr) = ~~ 2 (syst. dominated) , (B5) 

as would be expected intuitively. 

Let us divide now the spectra in n bins of width ATf. p = ATi} i= i t „^ n and p' = 
{p[, ATj}i = i r .. in , with J2i Pi^Ti = J2i Pi^Ti = 1. Then the shift in the height of the i-th bin 
associated to the linear deformation in Eq. (|B1|) is: 



Afi = d-n = fHp Ti W , (B6) 

where Tj is the average value of T in the z-th bin. 

If counting statistics dominates the errors, the fractional uncertainty <Ji/pi of the i-th 
bin height is Uj/pj = 1/^NpiATi, and the total x 2 °f the histogram differences, xlist = 
^j(Apj/(jj) 2 , is easily derived: 

* 2 

xLt = P 2 ^y 2 N , (B7) 

where a 2 = J2i PiATj(T, — (T)) 2 is just a discretized estimate of the variance a 2 , and therefore 
ex 2 — cr 2 . One gets the desired proof by comparing Eqs. ( |B4| ) and (|B~7|): 

xLt - x\t) ( stat - dominated) . (B8) 

If systematic errors dominate, the corresponding analysis of a binned spectrum is some- 
what trickier. Let the error be represented by an overall uncertainty op of the slope param- 
eter p. From Eq. (|B6|), this uncertainty propagates to an error cr^ = appiiTi — (T))/(T) of 
the i-th bin. The formula Xhist = J^ii^Pi/^i) 2 naively (and incorrectly) applied to this case 
would give Xhist = n P 2 1°% = n X(T)i with xfr) gi ven by Eq. ( |B5| ) . However, a systematic 
shift in the slope P produces completely correlated bin errors: corr(i, j) = 1. Therefore, 
only one out the n bin residuals is independent, and the "effective" number of bins to be 
considered in the x 2 is 1 0; so that Xhist = P 2 / a 1 an d : 



9 The reader more experienced in statistical analyses may have noticed that, in this case, the 
square error matrix including the error correlations would have rank 1 and not n, signaling that 
only one bin error is independent. 
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Xhisto 



(syst. dominated) 



(B9) 



Equations (B8) and (|B9| ) show that, if the errors are dominated by statistics or by a 
single systematic uncertainty, the use of the integrated variable (T) is as informative as a 
spectrum histogram, provided that the spectral deformations are linear in T. 
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TABLES 



TABLE I. The percentage la errors from different ingredients that affect the standard predic- 
tions, (T e ) = 7.658 MeV and Rqc/Rnc = 1.882. The numbers given are for N C c = 5000 events 
above threshold (T m i n = 5 MeV), and ecc = L £nc = 0.5. Uncertainties due to the backgrounds 
are neglected. The approximate cancellation of the correlation of total the errors may not be as 
strong for the actual SNO error budget. 



Error component a((T e )) % o{Rqq/ -Rnc) % Correlation 



Neutrino spectrum 


0.38 


0.43 


+1 


Cross section 


0.14 


0.53 


~ 


Statistics 


0.33 


3.09 





Energy resolution 


0.31 


0.47 


-1 


Energy scale 


0.68 


1.81 


+1 


Efficiency 


~ 


2.00 


~ 



TOTAL 0.91 4.18 0.32 



TABLE II. Deviations of (T e ) and -Rcc/-^NC from the standard model (electron flavor con- 
served) predictions. The results are shown for representative neutrino oscillation scenarios in units 
of standard deviations (cr). In the last two columns, the combined x 2 f° r the CC-shape and CC/NC 
ratio tests is calculated with and without the correlation of the total errors (p = 0.32), and the 
deviation is given as: Dev. (a) = \J~x^- Uncertainties due to the backgrounds are neglected. 

CC-shape test CC/NC test Combined tests 
(T e ) (MeV) Rcc/Rnc P = 0.22 p = 



Scenario 


Acronym 


Value 


Dev. (a) 


Value 


Dev. (a) 


Dev. (cr) 


Dev. (, 


Standard 


STD 


7.658 




1.882 








Small Mixing Angle (MSW) 


SMA 


7.875 


3.1 


0.639 


15.7 


17.9 


16.0 


Large Mixing Angle (MSW) 


LMA 


7.654 


0.0 


0.422 


18.5 


19.5 


18.5 


Vacuum Oscillations 


VAC 


8.361 


10.0 


0.411 


18.6 


25.1 


21.1 



TABLE III. The dependence of (T e ), Rqc/Rnc, and of their errors, on the threshold energy, 
T m i n . The correlation of errors is given in the last column. The results shown assume 5000 
CC events collected above threshold. Uncertainties due to the backgrounds are neglected. The 
acronyms for the scenarios (STD, SMA, LMA, VAC) are the same as in Table 0. 



T min Average energy (T e ) (MeV) 

(MeV) STD ± la SMA LMA VAC 



4.0 


7.234 ± 0.079 


7.533 


7.228 


8.101 


5.0 


7.658 ± 0.070 


7.875 


7.654 


8.361 


6.0 


8.187 ±0.063 


8.337 


8.184 


8.678 


7.0 


8.798 ± 0.058 


8.897 


8.796 


9.110 



CC/NC ratio, R CC /Rnc 




Correlation 


STD ± la 


SMA 


LMA 


VAC 


P 


2.177 ±0.088 


0.712 


0.489 


0.440 


0.23 


1.882 ± 0.079 


0.639 


0.422 


0.411 


0.32 


1.509 ± 0.067 


0.534 


0.338 


0.369 


0.43 


1.107 ±0.056 


0.409 


0.248 


0.307 


0.52 
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TABLE IV. The energy threshold (T m i n ) dependence of the deviations, £, of the predictions 
with neutrino oscillations from the standard predictions. The entries in the table, £((T e )) and 



?l-rCCC/-TLNCj 


, axe in 


units of standard deviations. 


The acronyms 


for the oscillat 


ion scenarios (SMA, 


T TV J A T 7" A /~1 \ 

LMA, VAC) 


are trie 


same as 


in Table ||. Uncertainties due to the backgrounds are ne£ 


dected. 


T ■ 

± mm 




f«r e » 




Z(Rcc/Rnc 


) 


Combined, £ = 




(MeV) 


SMA 


LMA 


VAC 


SMA 


LMA 


VAC 


SMA 


LMA 


VAC 


4.0 


3.8 


0.1 


11.0 


16.6 


19.2 


19.7 


18.3 


19.8 


25.4 


5.0 


3.1 


0.0 


10.0 


15.7 


18.5 


18.6 


17.5 


19.5 


25.1 


6.0 


2.4 


0.0 


7.8 


14.6 


17.4 


17.0 


17.4 


19.2 


23.9 


7.0 


1.7 


0.0 


5.4 


12.5 


15.3 


14.3 


15.8 


17.9 


20.7 



TABLE V. The CC event rates as a function of the CC energy threshold and the neutrino 



oscillation solution. 



T ■ 

± mm 


STD 


SMA 


LMA 


VAC 


(MeV) 


(SNU) 


(SNU) 


(SNU) 


(SNU) 


4 


6.9 


2.3 


1.6 


1.4 


5 


6.0 


2.0 


1.3 


1.3 


6 


4.8 


1.7 


1.1 


1.2 


7 


3.5 


1.3 


0.8 


1.0 
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FIGURES 



FIG. 1. (a) Total CC cross section as calculated by Kubodera and Nozawa (KN), Ying, Haxton 
and Henley (YHH), and Ellis and Bahcall (EB), slightly improved, (b) Total NC cross section as 
calculated by Kubodera and Nozawa (KN), and Ying, Haxton and Henley (YHH). 

FIG. 2. Comparison of the CC differential cross-section at various energies and scattering 



angles. Solid: Kubodera and Nozawa [27|. Dotted: Ellis and Bahcall [30|, slightly improved 



FIG. 3. The normalized electron spectrum, with and without inclusion of the energy resolution 
function. 

FIG. 4. Three standard deviation departures from the standard electron spectrum (solid line) 
due to neutrino-related and detector-related errors. 

FIG. 5. Neutrino Oscillation Scenarios: (a) survival probabilities for oscillation test cases; (b) 
effect of neutrino oscillations on neutrino spectrum at earth; (c) effect of neutrino oscillations on 
normalized electron spectrum at SNO. Area under curves = 1; (d) effect of neutrino oscillations 
on electron spectrum at SNO. Area under curves = Nqc/^nc- Labels: STD = standard (no 
oscillation); SMA = small mixing angle (MS W); LMA = large mixing angle (MSW); VAC = vacuum 
oscillation. See the text for details. 

FIG. 6. Ratios of the normalized neutrino spectra for different oscillation scenarios. The 
normalized spectra are displayed in Fig. ||c. Labels as in Fig. 5. 

FIG. 7. Values of the characteristic CC-shape variable, the average energy (T e ), and of the 
CC/NC ratio, Rcc/Rnc, together with 3a error bars. Uncertainties due to the backgrounds are 
neglected. Labels as in Fig. 5. 

FIG. 8. Iso-sigma contours (a = vx) f° r the combined CC-shape and CC/NC test, for the 
representative oscillation cases shown in Fig. 5 and discussed in the text. Uncertainties due to the 
backgrounds are neglected. For values of the iso-sigma distance ftf(a) 3, the number of standard 
deviations is only a formal characterization; the tail of the probability distribution is not expected 
to be Gaussian for very large values of N(cr). Labels as in Fig. 5. 

FIG. 9. A comparison of the 8 Li beta-decay spectrum and the standard electron spectrum 
from 8 B neutrino absorption, as a function of the electron kinetic energy above the standard SNO 
threshold (5 MeV). The spectra shown are both theoretical: the effects of finite energy resolution 
are not included. The area under the curves is normalized to unity. 
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FIG. 10. (a) Theoretical 8 Li spectrum and its 3a uncertainties, (b) Experimental determina- 
tions of the 8 Li spectrum, (c) Experimental data with an allowance for a linear recalibration of 
the energy. See the text for details. 
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FIG. 1. (a) Total CC cross section as calculated by Kubodera and Nozawa (KN), Ying, Haxton 
and Henley (YHH), and Ellis and Bahcall (EB), slightly improved, (b) Total NC cross section as 
calculated by Kubodera and Nozawa (KN), and Ying, Haxton and Henley (YHH). 
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FIG. 2. Comparison of the CC differential cross-section at various energies and scattering 



angles. Solid: Kubodera and Nozawa [27|. Dotted: Ellis and Bahcall [ 30 1 , slightly improved 
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FIG. 3. The normalized electron spectrum, with and without inclusion of the energy resolution 
function. 
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FIG. 4. Three standard deviation departures from the standard electron spectrum (solid line) 
due to neutrino-related and detector-related errors. 
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FIG. 5. Neutrino Oscillation Scenarios: (a) survival probabilities for oscillation test cases; (b) 
effect of neutrino oscillations on neutrino spectrum at earth; (c) effect of neutrino oscillations on 
normalized electron spectrum at SNO. Area under curves = 1; (d) effect of neutrino oscillations 
on electron spectrum at SNO. Area under curves = -/Vcc/^Vnc- Labels: STD = standard (no 
oscillation); SMA = small mixing angle (MSW); LMA = large mixing angle (MSW); VAC = vacuum 
oscillation. See the text for details. 
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FIG. 6. Ratios of the normalized neutrino spectra for different oscillation scenarios. The 
normalized spectra are displayed in Fig. ||c. Labels as in Fig. 5. 
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FIG. 7. Values of the characteristic CC-shape variable, the average energy (T e ), and of the 
CC/NC ratio, Rcc/Rnc, together with 3a error bars. Uncertainties due to the backgrounds are 
neglected. Labels as in Fig. 5. 
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FIG. 8. Iso-sigma contours (<r = y^?) f° r the combined CC-shape and CC/NC test, for the 
representative oscillation cases shown in Fig. 5 and discussed in the text. Uncertainties due to the 
backgrounds are neglected. For values of the iso-sigma distance N(a) » 3, the number of standard 
deviations is only a formal characterization; the tail of the probability distribution is not expected 
to be Gaussian for very large values of N{a). Labels as in Fig. 5. 
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FIG. 9. A comparison of the 8 Li beta-decay spectrum and the standard electron spectrum 
from 8 B neutrino absorption, as a function of the electron kinetic energy above the standard SNO 
threshold (5 MeV). The spectra shown are both theoretical: the effects of finite energy resolution 
are not included. The area under the curves is normalized to unity. 
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FIG. 10. (a) Theoretical 8 Li spectrum and its 3d uncertainties, (b) Experimental determina- 
tions of the 8 Li spectrum, (c) Experimental data with an allowance for a linear recalibration of 
the energy. See the text for details. 
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